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Abstract 

We study the influence of network topology and connectivity on the synchronization 
properties of chaotic logistic maps, interacting with random delay times. Four dif- 
ferent types of topologies are investigated: two regular (a ring-type and a ring-type 
with a central node) and two random (free-scale Barabasi-Albert and small-world 
Newman- Watts). The influence of the network connectivity is studied by varying the 
average number of links per node, while keeping constant the total input that each 
map receives from its neighbors. For weak coupling the array does not synchronize 
regardless the topology or connectivity of the network; however, for certain connec- 
tivity values there is enhanced coherence. For strong coupling the array synchronizes 
in the homogeneous steady-state, where the chaotic dynamics of the individual maps 
is suppressed. For both, weak and strong coupling, the array propensity for synchro- 
nization is largely independent of the network topology and depends mainly on the 
average number of links per node. 
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The emergence of dynami cal order in complex sy stems has been widely stud- 



ied during the last years (jPikOll ; lMan04t lBoc06[ ). One important issue is the 



propensity for synch ronizati on of networks of dynamical elements. In this con- 
text, coupled maps are excellent tools for understanding the mecha- 
nisms of emergency of synchrony and collective behavior in complex systems 
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composed of mutually coupled nonlinear units. Not only from an academi- 
cal point of view but also from an applied perspective, cooperative behavior 
arises in many fields of science and classical examples include the onset of 
rhythmic activity in the brain, the flashing on and off in unison of populations 
of fireflies, the emission of chirps by population s of cri ckets and the synchro- 



nization of laser arrays and Josephson junctions ([PikOlf ). Coupled map lattices 
have proven to be a useful tool because by simplifying the dynamics of the 
individual units it is possible to simulate large ensembles of coupled units. 

Network synchronizability and its relation with the to pology has recently re- 



ceived a great deal of attention. F. Atay et al. (JAtaOJ) found, in the case of 
fixed (constants) delays that scale-free and random networks exhibit better 
syn chroniza tion properties than regular networks. More recently, A. Motter et 
al. (|Mot05l ) identified that the synchronization of complex networks follows a 
diffusive mechanis m where the mean connectivity plays a key role. In a previ- 



ous paper we investigated the relation of the topology with the ability 



to synchronize under the presence of random delays. We found that in this 
case the synchronization properties depend largely on the mean connectivity 
of the network. However, the topology does not play an important role. The 
aim of this paper is to further investigate this point. 

The evolution equations for iV coupled logistic map with random delays are 

e N 

Xi(t + 1) = (1 - e)f[xi(t)] + //,,/ - Tij )]. (1) 

°i j= i 



Here t is a discrete time index, i is a discrete spatial index (i — 1 . . . N), f(x) = 
ax(l — x) is the logistic map, the matrix 77 = (77^) defines the connectivity 
of the array: rjij = rjji = 1 if there is a link between the ith and jth nodes, 
and zero otherwise, e is the coupling strength and is the delay time in the 
interaction between the ith and jth nodes (the delay times and need 
not be equal). The sum in Eq. (1) runs over the 6, nodes which are coupled 
to the zth node (bi = J2j Vij)- The normalized pre-factor 1/bi means that each 
map receives the same total input from its neighbors. 



In a previous work ( Mas05l ) we found that if the delays random enough, 



for adequate coupling strength the array synchronizes in the spatially homo- 
geneous stead y-state, X j(t) = x for all i, where x is the nontrivial fixed point, 
Xq — 1 — 1/a ( Mas05l ) . This synchronization behavior is in contrast with the 



behavior with fixed delays [if Ty — Tq V i, j, the array synchronizes in a 
spatially homogeneous time- depe ndent st ate, where the dynamics is either pe- 



riodic or chaotic depending on tq (JAtaOJ)], and can be understood in terms of 



the anal ogy between gl obally coupled maps and a single map with a external 
driving ( CosOl : MarOd ) . 
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Fig. 1. (Color online) Synchronization regions for the four different networks consid- 
ered. The density plots represent the parameter a 2 as a function of e and b (N = 500 
and a = 4.). 

To investigate the influence of the topology we consider four networks, two 
of them are regular and the other two are random. The regular ones are a 
ring of nearest-neighbor elements while in the second one we added a central 
node connected to all other nodes. The random networks consist of a scale free 
network constructed according to the Barabasi- Albert method and, concerning 
the last one, we use the small-world topology proposed by Newman and Watts. 
To characterize the transition to synchronization we use the indicator cr 2 = 
l/N(J2i[xi(t) — {x)] 2 ) t , where (.) denotes an average over the elements of the 
array and (.) t denotes an average over time. If the array synchronizes in a 
spatially homogeneous state, Xi(t) = Xj(t) V i, j, and, obviously, cr 2 = 0. 

We consider Gaussian distributed delays: r%j = r + near(c£), where c is a 
parameter that allows varying the width of the delay distribution; £ is Gaus- 
sian distributed with zero mean and standard deviation one; near denotes the 
nearest integer. Depending on r and c the distribution of delays has to be 
truncated to avoid negative delays. Since the focus of this paper is the influ- 
ence of the array topology and the connectivity we keep the random delays 
Gaussian distributed with (t^) ~ tq = 5 and c = 2. The numerical results are 
summarized in figures 1-3. 
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Fig. 2. a 2 as a function of the number of neighbors b for the weak coupling regime 
and different topologies: nearest-neighbors with central node (o), scale-free (□), 
nearest-neighbors (O) and small- world (A). Parameters are: N = 500, a = 4, (a) 
e = 0.02, (b) e = 0.04 and (c) e = 0.06 

In Fig. 1 we can see density plots of o 2 as a function of the mean number of 
links per node, b = l/NJ2iLi h and e. The four different panels correspond to 
the different networks mentioned above. Despite the differences for small num- 
ber of neighbors, we observe that the synchronizability is largely independent 
of the topology. Furthermore, the similarity between Figs. 1(c) and (d) clearly 
suggest that the synchronizability does not depend on the average path length 
(denned as the distance between two nodes, averaged over all pairs of nodes). 

For weak coupling (roughly speaking, e < 0.4), the array does not synchronize 
regardless the number of neighbors and topology. However, there are worth 
noting different behaviors depending on the value of e. To gain additional 
insight, it is shown in Fig. 2 the value of a 2 as a function of b for different values 
of e. We observe that for e = 0.02 (a) the value of a 2 decreases monotonously 
as b increases, to a limiting non-zero value. For e = 0.04 (b) there is still 
no synchronization for any connectivity value, but there is a non-monotonous 
dependence of a 2 with b, that reveals the existence of an optimal number 
of neighbors for which there is an enhancement of the array propensity for 
synchronization. Finally, in (c) we observe that, for e = 0.6, in sharp contrast 
to (a), a 2 increases monotonously before reaching a non-zero limiting value. 

On the other hand, for strong enough coupling (e > 0.4) the array synchronizes 
in the homogeneous steady state if the number of the neighbors is large enough. 
In figure 3 we observe a 2 as functions of b for different values of e. Depending 
on the value of e there is a minimum value of neighbors required to synchronize 
the array. Moreover, this critical value decreases with increasing e. 



In summary, we studied the synchronization of coupled maps in complex net- 
works with time-delayed interactions focusing on the influence of array con- 
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Fig. 3. a 2 as a function of the number of neighbors for the strong coupling regime 
and different topologies, (a) e = 0.6, (b) e = 0.8 and (c) e = 0.9. Same symbols and 
parameters values as Fig. 2. 

nectivity and topology. For weak coupling no synchronization was found, but 
an enhancement of the synchronization propensity was observed for particular 
connectivity values. On the other hand, for strong coupling there is synchro- 
nization provided that the number of neighbors was large enough. 
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